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ABSTRACT 

Aims. Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different 
types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models 
at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar 
metallicity, where the effects of rotation are accounted for in a homogeneous way. 

Methods. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z = 0.014, spanning a 
wide mass range from 0.8 to 120 Mq. For each of the stellar masses considered, electronic tables provide data for 400 stages along 
the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed 
data basis for comparisons with the observations. The rotating models start on the zero-age main sequence (ZAMS) with a rotation 
rate Vim/fcrit = 0.4. The evolution is computed until the end of the central carbon-burning phase, the early asymptotic giant branch 
(AGB) phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial 
abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar 
neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as 
stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass 
star models. 

Results. The present rotating models provide a good description of the average evolution of non-interacting stars. In particular, they 
reproduce the observed main-sequence width, the positions of the red giant and supergiant stars in the Hertzsprung-Russell (HR) 
diagram, the observed surface compositions and rotational velocities. Very interestingly, the enhancement of the mass loss during 
the red-supergiant stage, when the luminosity becomes supra-Eddington in some outer layers, help models above 15-20 Mq to 
lose a significant part of their hydrogen envelope and evolve back into the blue part of the HR diagram. This result has interesting 
consequences for the blue to red supergiant ratio, the minimum mass for stars to become Wolf-Rayet stars, and the maximum initial 
mass of stars that explode as type II-P supernovae. 
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1. Introduction least in the framework of the single star scenario. Changes in 

the abundances are observed at the surface of stars, which im- 

H ; Even if stellar physics has already a well-established history, pUes that mixing pr ocesses operate in radiative zones (see e.g. 

P. , which has allowed this area of research to reach a high level of iPrzybilla et al.|201Q). The origin of the valuation with metallicity 
sophistication, many surprises and c hallenges contmue to stim- in the numbe r ratio of blue to red supergiants still lacks a general 
ulate ongoing studies. For instance, lAsplund et all (l2005h pro- explanation (Mevlan & Maeder| ll982|; |Langer & Maeder||l995|; 
posed a revision of the present-day surface solar abundances Eggenberger et al. 2002). These have been only a few illustra- 
which restores the compatibility of the composition of our star tions the fact that, although many aspects of the evolution of 
with the one observed in B-type stars in its vicinity. However stars are well understood, important problems remain unsolved, 
these abundances pose some difficulties when attempting to ex- Significant progress is now achievable using observations, 
plain the heUoseismic observations, which may indicate that Multi-object specti'ographs can acquire data for lai-ge numbers of 
some pieces of physics a re still missing in solar models (e.g. stars, unveiling an unprecedentedly detailed view of how some 
iTurck-Chieze et al.1 120T1I and references therein). New mass- stellar features vary as a function of mass, age, and metalhcity. 
loss rates (in general lower than the previous ones) have been Techniques such as asteroseismology, spectropolarimetry or 
propos ed for the main-seq uence (M S) and Wolf-Rayet (WR) interferometry provide complementary insight, making it possi- 
phases dBouret et a l. 2005; Fullerton e t al. 2006) making it more ble to obtain information on features such as the size of the con- 
difficult to explain for instance the formation of WR stars, at vective cores, the internal rotation laws, the surface magnetic 



2 



Ekstrom et al.: Grids of stellar models with rotation 



fields, and the shapes of the stars and their winds, which were 
completely out of reach only a few years ago (ICunha et alJ2007t 
[Walder et al. 2011). 

In the near future, the ESA new satellite GAIA will provide 
data for large fractions of stars in our Galaxy to help enha nce our 
knowledge of its evolutionary history (see e.g. ICacciar 

Stellar models play a crucial role in extracting from 
observations the physical quantities needed to enhance our 
knowledge. Therefore, improvements of the observing fa- 
cilities should be accompanied by improvements in the 
stellar models. In the past decade, the effects of ax- 
ial rotation have been studied extensively by many groups 
dHeeer & Langejl2000l: iMeynet & Maedeill2000t iPalacios et al 



Table 1. Chemical initial abundances of the models given in 
log(X/H) +12 units and in mass fraction. 



Element Isotope 



2003 



iTalon & Charbonnef 20051; Denissenkov & P insonneault 
20071 Chieffi & Limongi, in prep.). It does appear that rotation 
has an impact on all outp uts of stellar models (Meynet & Maede^ 
I2OOOI: iMaeder & Meynetil2010l) . such as lifetimes, evolutionary 
scenarios, or nucleosynthesis, thus on many outputs of both 
population synthe sis models and ch emical evolution models 
of star clusters (Decressin et al. "2007), starbursts, and galaxies 
dChiappini et al. 2006; Vazquez et al. 2007). 

Grids including rotation have so far been focused on some 
limited mass intervals. Since the grids published by the different 
research groups differ in terms of their input physical ingredients 
(overshooting, mass loss, nuclear reaction rates,...) or the way in 
which the effects of rotation are taken into account (expressions 
of the diffusion coefficients, treatment of the angular momentum 
equation), their results cannot be combined to provide a coher- 
ent and homogeneous picture. In the present series of papers, our 
ambition is to construct a new database of evolutionary models 
incorporating the main improvements made in recent years to 
stellar ingredients and accounting in particular for the effects of 
axial rotation. We take care to compute the models in a homoge- 
neous way allowing us to provide a coherent and consistent view 
over an extended mass and metallicity range. In this first paper, 
we focus on solar metallicity models. 

The paper is structured as follows: the physical ingredients of 
the models are discussed in Sect. |2] The computed stellar mod- 
els and the electronic tables, which we make publicly available 
via a web page, are presented in Sect. [3] The properties of the 
non-rotating (or initially slowly rotating) tracks are discussed in 
Sect, m while the rotating tracks are presented in Sect. |5] We 
briefly conclude in Sect. |6l 

2. Physical ingredients of the models 

2.1. Abundances, opacities, and equation of state 

The initial abundances of H, He, and metals are set toX = 0.720, 
Y = 0.266, and Z - 0.014. We obtained these initial abundances 
by calibrating a 1 Mq model including atomic diffusion to repro- 
duce the abundances measured at the surface of the Sun, as well 
as the solar radius and luminosity after 4.5 7 Gyr. The mixture of 
heavy elements is assumed to be that of lAsplund e t al.l (^20051) 
except for the Ne abundance, which is taken from [Cunha et^ 
('2006'). Interestingly, the abundances correspond to the mea- 
sured abu ndances in young mas sive stars of the solar neighbour- 
hood (see lAsplund et al.ll2009l their Table 4). The detailed el- 
ement abundances are prese rited in Table [T] We note that the 
isotopic ratios are taken from lLoddersI (l2003l) . 

The opacities are generated with the OPAL tooQ (based 
on esias & Rogers, ,1996.) for this particular mixture. They 
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are complemented at low temperatures by the opacities from 
jierguson et al. (2005) adapted for the high Ne abundance. 

Solar-type models with M < 1 .25 Mp^ are com puted with 
the OPAL equation of state ("Rogers & Navfonovl200 2). For the 
higher mass models, the EOS is that of a mixture of perfect gas 
and radiatio n, and account for pa rtial ionisation in the outermost 
layers, as in lSchaller et alJ (Il992h . and for the partial degeneracy 
in the interior in the advanced stages. 

2.2. Nuclear reaction rates 

The nuclear reaction rates are generated with the NetGe n toofl 
The y are taken mainly from the Nacre database ( Angul oet al.l 
1999), although some have been redetermined more recently and 
updated. We list them below, with a comparison to NACRE val- 
ues and a short description of the effects on stellar evolution. 

The ra te of '"^N(p , y)'^0 is taken from 
Mukhamedzh anov et al.1 (120031) . It is about half the NACRE 
value for temperatures below 10^ K, and compares we l l with 
ot her determinations such as lAngulo & Descouvemonll (l200Th 
or iLemut et alJ (l2006h . This reaction is the slowest of the CNO 
cycle. In the low-mass domain, the effects of lowering t his rate 
have been studied by Imbriani et al. (2004) and We iss et al.l 
(l2005h . They describe a slower H-burning process, and shal- 
lower temperature profiles leading to a more extended and 
slightly hotter core. The turn-off point is shifted towards a 
hi gher luminosit y. In the interm ediate-rn ass domain, the studies 
of iHerwig et all (12006) and 1 Weiss et alJ (12005,) show that with 



' |http : //adg . llnl . gov/Research/OPAL/opal . html 
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lower rates, the MS evolution occurs at a higher luminosity, and 
that later, the blue loops during core He burning get significantly 
shorter. 

The rate of '*He(aa,y)'2c is taken from lFvnbo et alj(l2005h . 
This rate is around twice the NACRE value below 25 ■ 10^ K, 
sUghtly lower than NACRE between 25 ■ 10'' and 10** K, similar 
to NACRE between 10*^ and 10'' K, and a tenth of the NACRE 
value above lO** K. The effects of this new rate on stellar evolu- 
tion were studied bv lWeiss et"!!] (l2005h and lHerwig et al.1 (|2006|) 
for intermediate-mass stars. The differences are extremely small 
in this mass range. In the massive star domain, a lower rate in 
the temperature range of core He burning is expected to lead 
to a slightly larger core at the end of the evol ution, with a low 
C/O ratio, and slightly lower '^C yields (Eks trom et al.1 1201 Oh . 
The uncer tainties i n this rate, as well as the ^'C(a, y)'^0 rate, is 
shown bv iTur et~al . (2009, 2010) to strongly affect the produc- 
tion of the medium-weight nuclei and the the weak .s-process 
nuclei. 

The rate of ^^Cia,yy^O is taken from iKunz eTal] (|2002|) . 
It is around 0.6-0.8 times the NACRE value below 6 ■ 10** K, 
and around 1.1-1.4 times the NACRE value above this tempera- 
ture. Both ! Weaver & WooslevI (119931) and llmbrianiet al.l dioiJlb 
explored the effects of varying this rate. A higher rate in the He- 
burning temperature range leads to larger cores, lower '^C, and 
higher '^O yields. 

Th e rate of ^^Ne(Q;, n)^^Mg is taken from I Jaeger et"al1 
(1200 Ih . It is around 0.6 times the NACRE value below 6 ■ 10** 
K. While there is little effect on stellar evolution, th e ratio of 
the isotopes of Mg are modified (iKarakas et al.ll20()6h . and less 
^^Mg and ^''Mg is produced. This reaction is a neutron source for 
the ^-process, and The et al. (2007) studied the consequences of 
reducing its rate, which they found would lead to a significant 
reduction in the i-process efficiency during core He burning. 

We note that the expected impact of the updated nuclear re- 
action rates is small and dwarfed by the effect of other modifica- 
tions. 

The models of massive stars (M > 9Mq) were computed by 
incorporating the NeNa-MgAl cycle. For this sub-network, the 
NACRE rates were also used, except for tw o reactions that were 
updated: ^'Ne(p y)^^Na (llliadis et al.' '2001) and ^^Ne{p, yf^Na 
dHale et al.ll2002l) . The updated reaction rates do not have a sig- 
nificant effect on the structure and evolution of the models. The 
effects of rotation and these new rates on the nucleosynthesis, in 
particular that of the i-process, will be discussed in a forthcom- 
ing paper (Frischknecht et al, in prep.). 

The energy loss in plasma, pair, and photo-neutrinos p ro- 
cesses are taken from lltoh etaP (Il989h and lltoh et all (Il996l) . 

2.3. Convection and diffusion 

The convective zones are determined with the Schwarzschild 
criterion. For the H- and He-burning phases, the convective 
core is extended with an overshoot parameter dova/Hp = 0.10 
from 1.7 Mq and above, 0.05 between 1.25 and 1.5 Mq, and 
below. If doYN exceeds the dimension of the convective core 
Rcc, then the total extension of the convective core is given 
by /?cc(l + doyei-/Hp). This procedure avoids any core-extension 
amplitude larger than the radius of the initial core. The value 
of the overshoot parameter was calibrated in the mass domain 
1.35-9 Mq to ensure that the rotating models closely reproduce 
the observed width of the MS band. 

In Fig.[T] we compare the MS bandwi dth of the present mod - 
els with that in the models published in ISchaller et alj (Il992h . 
Our new models concur to a slightly narrower MS width than 
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Fig. 1. Comparison of the MS band between the present rotat- 
ing models (solid b lack line) and the non-rotating models of 
ISchaller et al.l (Il992i) (long-short dashed blue line). 



Schaller et al.l (1199 2*). which provided a good fit to the observed 
MS width (see the discussion in Schaller et al. 1992). However, 
since stars have a distribution of initial rotation velocities, stars 
at the end of the MS phase will scatter around the limit shown 
by our moderately rotating models, more rapid rotators lying 
beyond it and slightly enlarging the MS width. A more de- 
tailed discussion and a comparison of the MS width obtained 
from our models with observations are provided in Sect. 15.11 
We note for now that the combined effects of rotational mixing 
and of an overshoot of 0. 1 mimic the effects obtained in mod- 
els with no rotational mixing and a stronger overshoot of 0.2 
(used in the 92's models). This is of course expected because ro- 
tational mixing also contributes to making the convective cores 
larger at a given evolutionary stage. This effect was discussed in 
Talon etal. (1997). 

The outer convective zone is treated according to the mixing 
length theory, with a solar calibrated value for the mixing-length 
parameter of the low-mass stars (aMLT = f^/Hp - 1.6467). For 
more massive stars, i.e. for stars with M > 1.25 Mq, the differ- 
ence in the EOS implies a slightly lower value for this parameter: 
ffMLT - 1-6. In Fig.|2] we see that the present models closely re- 
produce the positions of both the red giant branch and the red 
supergiants, supporting the choices made for the value of the 
mixing length. 

For the most luminous models, the turbulence pressure and 
acoustic flux need to be inc l uded in the treatment of the en- 
velope. As in ISchaller et al.' ('1992'), this is done according to 
Maeder (1987, see also Maeder 2009, Sect. 5.5), using a mix- 
ing length taken on the density scale: (Tmlt - (/Hp - {{a - 
6V)/Hp - 1. The use of (/Hp instead of (/Hp removes the 
risk of having an unphysic al density inversion in the enve- 
lope (IStothers & Chiniri973h . The side effect of this treatment 
is that the redward extension of the tracks in the Hertzsprung- 
Russell (HR) diagram i s reduced by 0.1-0.2 dex in Teff (see 
lMaeder&Mevnellll987h . which is why we apply it only to the 
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Fig. 2. Evolutionary tracks for rotating models in the red part 
of the HR diagram. In the low mass range, a few non-rotating 
tracks (dashed lines) are shown. The grey shaded areas indi- 
cate the observations in clusters and associations as given in 
lMaeder&Mevnet l('1989^ as well as the position of Galactic red 
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supergiants obtained bv .Levesque et al] ([2005t red circles). 



models with M > 40 Mq, which do not extend to the extreme 
red part of the HR diagram. 

Low-mass stellar models with M < 1 .25 Mq are computed 
with the effect of atomic dif fusion caused by both co ncentration 
and thermal gradients (see [Eggenberger et al.ll2008l for more 
details). 

2.4. Rotation 



The 
of 



treatment of rotation was developed in 
pers published previous ly by our group 



a series 

ot papers published previous ly by our group ([Maedeil 

19991 [19 97; Maeder & Mevnet l2000t iMaeder & ZahnI Il998t 



Mevnet & Maeder i1997l |2000). The interested reader may re- 
fer to these papers to get the full developments. We recall here 
the main aspects of this treatment. 

2.4.1. Diffusion coefficients, meridional circulation 

We first apply the shellular-rotation hypothesis, which postulates 
that in differentially rotating stars, the horizontal turbulence, i.e. 
the turbulence along an isobar, is very strong. This is expected 
because there is no restoring force in that direction, as the buoy- 
ancy force (the re storing force of the density gradient) acts in the 
vertical direction. IZahnl(ll992h relates the diffusion coeflicient to 



the viscosity caused by horizontal turbulence 



1 



£>h*Vh = - r \2Vir)-aUir)\, 



(1) 



where Ch is a constant of the order of 1 (here taken =1), V(r) is 
the horizontal component of the meridional circulation velocity, 
U{r) its vertical component (see below), and in this expression 

^ ^ 1 dln(/-^n) 

According to the von Zeipel theorem dvon Zeipe]|ll924t) . a 
rotating star cannot be locally in both hydrostatic and radia- 
tive equilibrium at the same time. This drives a large-scale cir- 
culation, known as the meridional c i rculation (ox Eddington- 
Sweet circulation, lEddingtonl IT925I: ISweei IT950l) . Its verti- 
cal component is to first order u(r, 0) - U{r) f 2(cos 0), where 
Piix) is the second Legendre polynomial The formulation of 
the radial ampli t ude U (r) was determined by IZahnI ((T992b and 
IMaeder &Z^ (1 19981) to be 

p 1 

U(r)^ 



L " 5 dtj 



(2) 



where Cp is the specific heat at constant pressure, V^d = 

is the adiabatic gradient, - m{\ - 5^")' ® - PiP is 
the ratio of the variation of the density to the average density on 
an equipotential. Both (/? and 5 arise from the generic equation of 
state y - and E*^ and /s^ are terms that depend 

on th e Q- and //-distributions respectively (see iMaeder & ZahnI 
for details on these expressions). 
We can then define an effective diffusion coefficient, D^s, 
which combines the effects of the horizontal diffusion and that 
of the meridional circulation as ( Chabover & Zahni,1992) 



eff ■ 



1 \rU{rr 
30 Dh 



(3) 



Differential rotation induces shear turbulence at the interface 
of layers that have different rotational velocities. A layer remains 
stable if the excess kinetic energy due to the differential rotation 
does not overcome the buoyancy force (known as the Richardson 
criterion). Shear arises when the thermal dissipation reduces the 
buoyancy force. The coefficient of diffusion by shear turbulence 
is determined bv lMaedeii([T997l) to be 



K 



shear 



f + (V,d - V,ad) 

Hp 



,9n dlnQ, , , 



(4) 



where K - 



4acT^ 



2i^p2(^^ is the thermal diffusivity, /energ is the fraction 
of the excess energy in the shear that contributes to mixing (here 
taken =1), and (V - V) is the difference between the internal 
non-adiabatic gradient and the local gradient. We note that (V - 
V) can be neglected in most cases. 



^ The Legendre polynomials are a set of orthogonal functions used to 
solve Laplace's equation V-O = in spherical coordinates when there 
is azimuthal symmetry dip = 0. The second-order Legendre polynomial 
is F2(x) = \ (3x2 _ ^ 
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2.4.2. Transport mechanisms 

The transport of angular mo mentum insid e a star is implemented 
following the prescription of Zahn'('1992"). This prescription was 
compl emented by Talon & Zahn (1997) and Ma eder & ZahnI 
(119981) . In the radial direction, it obeys the equation 

The first term on the right hand side of this equation is the di- 
vergence of the advected flux of angular momentum, while the 
second term is the divergence of the diffused flux. The coefficient 
D is the total diffusion coefficient in the vertical direction, tak- 
ing into account the various instabilities that transport angular 
mo mentum (convection, shea rs) . 

iChabover & (IT991 show that the horizontal turbu- 
lence competes efficiently with the advective term of the merid- 
ional circulation for transporting the chemical species. The hor- 
izontal flow tends to homogenise the layer in such a way that 
the transport of chemical species by both meridional circulation 
and horizontal turbulence can be computed as a diffusive pro- 
cess with the coefficient Defr calculated in Eq. [3] The change in 
the abundance for a given chemical element / in the shell with 
Lagrangian coordinate r is thus: 



(6) 



nucl 



where D is the same as in Eq. |5] The last term accounts for the 
change in abundance produced by nuclear reactions. 

2.4.3. Angular momentum conservation 

The final angular momentum content of the star (particularly 
the core content) is a key quantity for determinin g its final state 
(supernova, hypernova , gamma-ray burst, see e.g. lNomoto et al.l 
l2005HYoon etal.l2006h . The conservation of angular momentum 
during the whole stellar evolution is therefore of prime impor- 
tance. In the Geneva code, the angular momentum conservation 
is checked throughout the evolution of the star. During the com- 
putation of the models, there are two sources of variations in the 
total amount of angular momentum: 

- A numerical variation, owing to unavoidable inaccuracies 
in the resolution of the advection-diffusion equation of the 
angul ar momentum transport a nd the structure of the code 
itself (iKippenhahn et al.lll967l) : the envelop^ is "floating" 
over the interior layers, and is assumed to rotate at the same 
angular velocity than the most superficial shell of the inte- 
rior. Its angular momentum content is therefore imposed by 
what happens in the interior. For the total angular momen- 
tum of the star to remain constant, a correction needs to be 
applied (see below). 

- An angular momentum loss due to the mass loss of the rotat- 
ing star, where the stellar winds carry away some amount of 
momentum. 

Knowing the mass-loss rate M of the star (as a func- 
tion of stellar parameters such as the luminosity, the effective 
temperature, and the chemical composition of the surface, cf. 



In our models, the envelope is the region above a given mass co- 
ordinate in which the luminosity is considered as constant, and where 
partial ionisation is accounted for. 



Sect. 12.6.11 ). it is possible to compute the amount of angular mo- 
mentum AXwinds carried away by the stellar winds during the 
current time step At 



A-Cwinds = -MQs^Af, 



(7) 



where £ls is the angular velocity of the surface, and r, the stel- 
lar radius. We also assumed that the mass loss is spherically 
symmetric, even if some anisotropy may develop. iGeorgv et al.l 
(1201 lb demonstrated that this assumption leads to only small er- 
rors, particularly for stars that remain far from the critical rota- 
tioij^ as in this work. Starting from the total angular momentum 
X-ini of the model at the time step n - 1, the angular momentum 
conservation ensures that the expected angular momentum at the 
end of the n* time step should be 



fin 



-^ini A^. 



winds - 



(8) 



As discussed above, the final angular momentum obtained 
after the whole computation of the structure differs usually from 
the expected one -C^^- We thus need to correct the angular mo- 
mentum of the model by an amount AXcorr - -C^^ - -C^^''- The 
correction for the angular velocity of the /* shell is of the form 
of fi™" - Q.f(l +qcon) where Q°'' is the angular- velocity of 
the J* shell obtained by the numerical computation, and ^con the 
coiTection factor, given by 



(9) 



where Xe and X.i are the angular momenta of the envelope and 
the shell, re spectively. A d etailed description of the method 
can be found in lGeorgvl(l20Toh . 

The choice of A^corr is fixed, and corresponds roughly to the 
zone that can be reached by the advection-diffusion of the angu- 
lar momentum during one characteristic time step. 



2.5. Magnetic field 

We accounted for neither a dynamo mechanism in the stellar in- 
terior, nor any strong fossil field that would impose solid body 
rotation. Our present knowledge of the dynamo theory is too un- 
certain (jZahn et al. 2007), and the observational constraints from 
spectropolarimetry remain too weak t o account for these effects 
in a reliable way. According to .Maeder & Mevnetl (l2005h . an in- 
ternal magnetic field produces a strong internal coupling, which 
keeps the star at a higher surface rotational velocity throughout 
the whole MS. This may provide an interesting test of whether a 
significant internal magnetic coupling is present. 

We note that the main reasons for introducing the effects of 
internal magnetic fields that would enforce solid body rotation 
during the MS phase come from essentially two directions: to 
provide first an explana tion for the flat rotational profile in the 
Sun (Brown et al. 1989 : iKosovichev et alJll997t ICouvidat et all 
I2003i: lEggenberger et all l2005h . and second, a more efficient 
mechanism for slowing down the cores to explain th e observed 
rotation rates of young pul sars (Heger et al.l l2005h and white 
dwarfs (WD, Suiis et al.ll2008) . The problems of the rotational 



' The critical velocity is reached when the gravitational acceleration 
is exactly counterbalanced by the centrifugal force. In the framework of 

the Roche model, one has Ucrit = ^| where /?pb is the polar radius 

at the critical limit. 
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profile of the Sun and the spin rate of pulsars might not be linked, 
since the physics of solar-type stars have peculiarities that are 
not shared by the massive ones. We note that in both ques- 
tions (solar profile and pulsar rotation), alternative models have 
been suggested to reconcile the models with the observations: 
the flat solar rotational profile co uld be due to angu l ar momen- 
tum transport by gravity waves ( Schatzman 1993"; IZahn et aP 
[1997; Talon et al. 2002; Talon & Charbonnel 2005), the young 
pulsars may have slower rotational velocities than predicted by 
stellar evolution models because of some efficient braking mech- 
anism occurring at the time of the supe rnova explosion (for 
instance, the propeller mechanism, see Illar ionov & SunvaevI 
ll975HAlparetalJl200lh or during the evolution of the neutron 
stars in the period preceding the time when the pulsa r is ob- 
served ( [Thompson et aLllIooi; IWooslev & H"ege^l200l . While 
kicks are not hkely to occur in WD formation, surface mag- 
netic braking (see next paragraph) might provide a solution to 
the problem of the lo w-mass WD spin rates (see the 1 Mq model 
of ISuiis et al.ll2()()8h . For all these reasons, we see no real com- 
pelling reasons to account for either a dynamo mechanism in 
the internal radiative zones or the effects of a strong fossil in- 
ternal magnetic field. In support of this view are the few but 
steadily increasing pieces of observational evidence now coming 
from asteroseismology for differential rotation in B-type stars 
(ISuarez et al ]l2007l: lThoulll2009t lKawaledl2()()9l) . 

Surface magnetic fields are however expected to produce a 
magnetic braking of stars that have a significant outer convective 
zone on the main sequence, hence magnetic braking has thus far 
been appli ed to models with M „ < 1 -7 Mq. We adopt the brak- 



ing law of iKrishnamurthi et al 



1997h 

braking constant to ensure that the 1 Mq rotating model repro- 
duces the solar surface rotational velocity after 4.57 Gyr. We did 
not account for surface magnetic braking in more massive stars, 
alth ough there is, in at least one case, evidence of magnetic brak- 
ing dTownsend et al.r20 10). Again, this effect needs to be studied 
in greater detail before being incorporated into the regular mod- 
els and we defer a more detailed study of the consequences of 
this effect to a future paper. 



2.6. Mass loss 

Mass loss is a key ingredient governing the evolution of stars. 
Including it in stellar models relies on prescriptions proposed by 
both observers and theorists and based on fundamental stellar 
parameters. The prescriptions implemented in our models are 
described below. 



2.6.1. Radiative mass loss 

On the MS, stars with masses below 7 Mq are computed at con- 
stant mass. Above 7_Mq, the radiative mass-loss rate adopted 
is from Vinket al. (2001). In the domains not covered by this 
prescription, we use de Jager etaL (1988). 

For red (super)giants, we use theH eimer sl (ll975Lll977h for- 
mula ( with T] - 0.5) for stars up to 12 Mq. The de Jager et alJ 
(Il988l) prescription is applied for masses of 15 Mq and above, 
to models with log(reff) > 3.7. For log(reff) < 3 . 7, we 
per form a hnear fit to th e data from ISvlv ester et al.l (Il998l) 
and Ivan Loon et al.) (Il999l) (see IC rowther 200(J. The WR stars 
are computed with Nugis & La mersI (12000 ) prescription, or the 
[Orafener & Hamann (2008) recipe in the small validity domain 
of this prescription. In som e cases, the WR mass-l oss rate of 
iGrafener & Hamamil (l2008h is lower than the rate of IVink et al.l 



(I2001h . In these cases, we use the IVink et al ] (I2001h prescrip- 
tion instead of that of Grafener & Hamann (2008). Both the 
Nugis & Lamers (2000) and Grafener & Hamann (2008 ) mass- 
loss rates account for some clumping effects (IMuiires et al.l 
i201 1 ) and are a fact or of two to three low er than the rates used 
in the normal case of lSchaller et al.l (Il992h grid. 

For rotating models, we applied to the radiative mass-loss 
rate the correction factor described in iMaeder & Mevnetl (|2000|) 



M(Q) = Fn ■ M(Q = 0) = ■ M,^d 

(i-r)i-' 



with Fq = 



(10) 



where F - L/Lem - kL/{47tcGM) is the Eddington factor (with 
K is the electron-scattering opacity), and a the force multiplier 
parameter depending on Teff. 

2.6.2. Supra-Eddington mass loss 

For some stellar models, particularly for massive stars (> 
15 Mq) in the red supergiant phase, some of the most external 
layers of the stellar envelope might exceed the Eddington lumi- 
nosity of the star LEdd = 4ncGM/K. This is due to the opacity 
peak produced by the variation in the ionisation level of hydro- 
gen beneath the surface of the star. The high opacity decreases 
the Eddington luminosity in these layers, possibly to fainter lev- 
els than the actual s tellar luminos ity, a situation that may have 
many consequences (lMaeded2009 )." To account for this unstable 
situation (which is not solvable with our hydrostatic approach), 
we artificially increase the mass-loss rate of the star (computed 
according to the prescription described in Sect. l2.6.TT i by a factor 
of 3, whenever the luminosity of any of the layers of the enve- 
lope is higher than five times the Eddington luminosity. 

Very interestingly, the time-averaged mass-loss rate of the 
20 Mq model in the red supergiant phase (i.e. when log(reif) is 
inferior to 3.65) is between log(M) = -4.8 and -4.6 [MQyr"'] 
(similar values are obtained for the rotating and the non-rotating 
model). These values are about one order of magnitude higher 
than th e mass-loss rates used for the sam e phase in Schaller et aT] 
(Il992h taken from lde Jager et all (119880 and are compatible with 
the more recen t mass- loss rates for red supergiants obtained by 
Ivan Loon et al. These new higher mass-loss rates dur- 

ing the red supergiant stage have important consequences for the 
blue to red supergiant ratio, the minimum mass of stars evolving 
into the Wolf-Rayet phase, and the type of supernova that these 
stars will produce. We briefly discuss these points in Sect. 15.51 

2.6.3. Mechanical mass loss 

We call mechanical mass loss, Mmech, the mass per unit time 
lost equatorially when the surface velocity of the star reaches 
the critical velocity, where Mmech is determined by the angular 
momentum that needs to be lost to ensure that the surface ve- 
locity remains subcritical. The matter might be launched into 
an equatorial Keplerian disk, a subject that will be studied in a 
forthcoming paper (Georgy et al. in prep.). 

The computational method used to estim ate the equatorial, 
mechanical mass loss is detailed in iGeorgvl (|2010). Since the 
initial velocity we have chosen for this grid is moderate, none of 
our models reach the critical limit during the MS. Only the 9 Mq 
model reaches it briefly during its Cepheid blue loop, losing at 
most 10"^ Mq in the process. In Georgy et al. (in prep), we will 
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consider the case of more rapid rotators, which reach the critical 
limit on the MS. 



3. The stellar models and electronic tables 

Models of 0.8, 0.9, 1, 1.1, 1.25, 1.35, 1.5, 1.7, 2, 2.5, 3, 4, 5, 
7, 9, 12, 15, 20, 25, 32, 40, 60, 85, and 120 Mq are presented. 
For each mass, we computed both a rotating and a non-rotating 
model. 

The rotating models start on the zero-age main-sequence 
(ZAMS) with a value of Vmilvcm - 0.4. We decided to com- 
pute the rotating models at a given Vini/L'clit ratio rather than a 
given initial equatorial velocity because the effects of rotation 
are linked to this ratio, and the critical velocity v„it depends on 
the mass of the model. For such a widespread mass domain, a 
given velocity would thus correspond to very different v/vait 
ratios between the different masses considered. It would even 
correspond to an over-critical velocity for some low-mass mod- 
els. The choice of Vini/vcrit - 0.40 is b ased on the peak of the 
velocity distribution of young B stars in lHuang et al.l (1201 Ol see 
their Fig. 6). This initial rotation rate corresponds to mean MS 
velocities of between 110 and 220 km s"' for the stars that are 
not magnetically braked (Mini > 1-7 Mq). These v alues are well 
within the range of those o bserved in the Galaxy dPufton et al.l 
l2006HHuang & Giesll2006h . 

The models are evolved up to the end of the core carbon 
burning (Mini ^ 12 Mq), the early asymptotic giant branch 
(2.5 Mq < Mini < 9 Mq), or the helium flash (Mini < 2 Mq). 

Electronic tables of the evolutionary sequences are available 
on the wet|3- For each model, the evolutionary track is described 
by 400 selected data points, each point corresponding to the state 
of the considered star at a given age. The points of different evo- 
lutionary tracks with the same line number correspond to simi- 
lar stages to facilitate the interpolation between the tracks. The 
points are numbered as follows: 

1: ZAMS 

2-84: H burning (first part) 

85: minimum of Teff on the MS 

86-109: overall contraction phase before the end of the MS 

110: Turn-off 

1 1 1-189: HR diagram crossing and/or pre-He-b core contraction 




Fig. 3. Isochrones from our rotating models. Log ages of 6.5 
(top) to 10.1 (bottom) are given, with steps of 0.1 dex. 



190: beginning of He burning 

191-209: He burning (first part) 

210-350: blue loop (if any, maximal extension on point 280) 
351-369: He burning (second part) 
370: core He exhaustion 
371-399: Churning 
400: last model 

For the stars that do not have a certain characteristic, we 
search for the closest (in mass) stellar model with that feature. 
The corresponding data point is then determined by the central 
abundance of the main fuel at that phase in this closest model. 
For example, the 12 Mq model does not evolve through the blue 
loop of a Cepheid. In this model, the point 280 is determined by 
the central mass fraction of He at the maximal extension of the 
loop in the 9 Mq model, which is the closest model with a loop. 
For the stars that do not go through all the burning phases, the 
last line of their evolution is repeated up to line 400, so we keep 
the same size for all the electronic tables. 



Isochrones are also provided (see Fig. O on the same web- 
page. We present the surface properties of the models at the 
given age, such as L, Teff, R, \ogg, surface abundances, rota- 
tional characteristics, as well as colour-magnitude values, such 
as Mboi, My, U — B, B — V, and B2 — V\ (co nversions from 
B6 hm-Vitenselll98ll: iFloweill 1 9 771: ISekiguchi & Fukugita 2005 
Ma lagninietal.l l 1986: Schmidt- Kaieilll982h . 

Table |2]presents the general characteristics of all the models 
at the end of each burning phase. After the initial mass, initial 
velocity, and mean velocity on the MS (col. 1 to 3), we give for 
each burning phase its duration (col. 4, 10, and 16), the total 
mass (col. 5, 11, and 17), the equatorial velocity (col. 6, 12, and 

18) , the surface He abundance in mass fraction (col. 7, 13, and 

19) , and the surface abundances ratios N/C (col. 8, 14, and 20) 
and N/O (col. 9, 15, and 21). 

When calculating the lifetimes in the central burning stages, 
we consider the start of the stage as the time when 0.003 in mass 
fraction of the main burning fuel is burnt. We consider that a 
burning stage is finished when the main fuel mass fraction drops 
below 10 ^ 

Table [3] offers a piecewise linear fit to calculate the MS du- 
ration from the mass of the star in the form log(rH) = A ■ 
\og(M I Mq) + B. The coefficients depend on the mass domain 
considered. As expected, the lifetime of stars decreases sharply 
with increasing stellar mass (A^ -3) in the low mass star regime 
and the dependence becomes milder (A^ - 1) in the massive star 
regime. The dispersion around the coefficients is also given in 
Table[3] and amounts at most to 0.025. 



4. Properties of the non-rotating models 

The present non-rotating tracks can be used for three different 
purposes: first they can be considered as describing the evolu- 
tion of slowly rotating stars, i.e. of stars rotating slowly enough 
for their outputs to be little affected by rotation. Second they al- 
^ See the webpage http : //obswww. unige . ch/Recherche/evol/-Ddt)SbS§eiparisons with non-rotating evolutionary tracks that were 
or the CDS database at http : //vizier . u- strasbg . f r/viz-bin/Vizi^Hita i shed previously. Third, since they were computed with ex- 
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log(N/H)+12 




4.5 4 3.5 

log (TJ [K] 

Fig. 4. Hertzsprung-Russell diagram for the non-rotating models. The colour scale indicates the surface number abundance of 
nitrogen on a log scale where the abundance of hydrogen is 12. Once the star has become a W NE type Wolf-Rayet st ar, the tracks 
are drawn with black dotted lines. The grey shaded area represents the Cepheid instability strip dTammann et alj|2003h . 



Ekstrom et al.: Grids of stellar models with rotation 



9 



actly the same physics as the rotating tracks, they offer a com- 
parison basis for studying the effects of rotation. 



Table 2. Properties of the stellar models at the end of the H-, He-, and C-burning phases. 



M 


Vini 






End of H-burning 










End of He-burning 










End of C-burning 












tu 


M 


t'surf 


Ysurf 


N/C 


N/O 


fee 


M 


i^surf 


Ysurf 


N/C 


N/O 


tc 


M 


fsurf 


Ysurf 


N/C 


N/O 


Mq 


km s 


-1 


Myr 


Mq 


km s"' 




mass fract. 




Myr 


Mq 


km s~' 




mass fract. 




kyr 


Mq 


km s~' 




mass fract. 




120 








2.672 


63.625 


- 


0.783 


94.89 


77.30 


0.308 


30.974 


- 


0.241 


0.00 


0.00 


0.008 


30.911 


- 


0.238 


0.00 


0.00 




389 


109 


3.178 


34.611 


1 


0.974 


70.15 


92.07 


0.350 


19.045 


3 


0.264 


0.00 


0.00 


0.030 


19.044 





0.264 


0.00 


0.00 


85 








3.024 


49.232 


- 


0.625 


113.01 


69.93 


0.348 


18.711 


- 


0.250 


0.00 


0.00 


0.031 


18.648 


- 


0.247 


0.00 


0.00 




368 


124 


3.715 


49.272 


4 


0.926 


81.20 


84.65 


0.319 


26.394 


19 


0.264 


0.00 


0.00 


0.013 


26.393 


23 


0.264 


0.00 


0.00 


60 








3.530 


36.261 


- 


0.493 


145.35 


60.81 


0.397 


12.565 


- 


0.279 


0.00 


0.00 


0.115 


12.495 


- 


0.273 


0.00 


0.00 




346 


138 


4.465 


38.536 


5 


0.820 


78.76 


49.35 


0.355 


17.981 


48 


0.277 


0.00 


0.00 


0.037 


17.981 


22 


0.277 


0.00 


0.00 


40 








4.439 


36.509 


- 


0.266 


0.29 


0.12 


0.478 


13.020 


- 


0.986 


63.68 


87.48 


0.066 


12.821 


- 


0.986 


33.14 


28.77 




314 


155 


5.698 


32.005 


20 


0.554 


12.73 


4.31 


0.424 


12.334 


72 


0.278 


0.00 


0.00 


0.123 


12.332 


55 


0.278 


0.00 


0.00 


32 








5.207 


30.107 


- 


0.266 


0.29 


0.12 


0.543 


11.151 


- 


0.647 


132.66 


41.13 


0.186 


10.922 


- 


0.986 


56.82 


77.97 




306 


187 


6.639 


28.195 


13 


0.441 


5.84 


1.81 


0.521 


10.127 


46 


0.280 


0.00 


0.01 


0.213 


10.125 





0.280 


0.00 


0.00 


25 








6.310 


24.180 


— 


0.266 


0.29 


0.12 


0.692 


8.757 


— 


0.537 


196.18 


6.69 


0.433 


8.289 


— 


0.831 


119.91 


72.30 




295 


209 


7.902 


23.604 


81 


0.341 


3.13 


0.86 


0.618 


9.700 





0.917 


102.03 


23.61 


0.269 


9.690 





0.926 


100.50 


25.12 


20 








7.740 


19.672 


- 


0.266 


0.29 


0.12 


0.874 


9.024 


- 


0.503 


59.43 


3.93 


1.219 


8.635 


- 


0.511 


81.52 


4.39 




274 


217 


9.506 


19.497 


176 


0.300 


2.42 


0.60 


0.846 


7.191 





0.737 


127.94 


14.62 


0.866 


7.179 





0.748 


127.30 


15.88 


15 








11.015 


14.812 


- 


0.266 


0.29 


0.12 


1.315 


13.340 


- 


0.307 


2.12 


0.53 


4.786 


13.249 


- 


0.339 


3.07 


0.74 




271 


201 


13.447 


14.701 


160 


0.289 


2.34 


0.53 


1.451 


11.079 





0.427 


9.40 


1.33 


1.374 


11.071 





0.429 


9.52 


1.35 


12 








15.330 


11.940 


- 


0.266 


0.29 


0.12 


2.124 


11.364 


- 


0.301 


1.84 


0.49 


7.927 


11.308 


- 


0.307 


1.98 


0.53 




262 


198 


18.369 


11.907 


192 


0.279 


1.89 


0.45 


2.098 


10.234 





0.349 


6.22 


0.91 


3.567 


10.224 





0.350 


6.27 


0.92 


9 








26.261 


8.993 


- 


0.266 


0.29 


0.12 


3.492 


8.802 


- 


0.282 


1.61 


0.42 


- 


- 


- 


- 


- 


- 




248 


188 


31.211 


8.988 


198 


0.274 


1.40 


0.37 


3.762 


8.521 


4 


0.337 


6.09 


0.86 


3.143 


8.517 


3 


0.337 


6.09 


0.86 


7 








41.721 


6.999 


- 


0.266 


0.29 


0.12 


6.919 


6.915 


- 


0.279 


1.52 


0.40 
















235 


178 


50.989 


6.999 


184 


0.270 


0.97 


0.29 


6.903 


6.892 


2 


0.322 


4.68 


0.76 














5 








88.193 


5.000 




0.266 


0.29 


0.12 


19.380 


4.959 




0.279 


1.48 


0.39 
















219 


166 


109.207 


5.000 


164 


0.268 


0.64 


0.21 


17.501 


4.946 


3 


0.318 


3.79 


0.69 














4 








152.082 


4.000 




0.266 


0.29 


0.12 


37.750 


3.975 




0.282 


1.50 


0.39 
















197 


159 


189.417 


4.000 


153 


0.267 


0.52 


0.18 


33.491 


3.967 


4 


0.319 


3.38 


0.66 














3 








320.585 


3.000 




0.266 


0.29 


0.12 


117.109 


2.988 




0.287 


1.53 


0.39 
















195 


147 


405.366 


3.000 


136 


0.267 


0.43 


0.16 


95.566 


2.984 


5 


0.323 


3.11 


0.63 














2.5 








537.935 


2.500 




0.266 


0.29 


0.12 


236.759 


2.487 




0.285 


1.44 


0.37 
















187 


141 


673.098 


2.500 


128 


0.267 


0.40 


0.15 


181.461 


2.486 


5 


0.322 


2.92 


0.59 














2 



186 



137 


1008.831 
1289.272 


2.000 
2.000 


124 


0.266 
0.267 


0.29 
0.38 


0.12 
0.14 


























1.7 



176 




133 


1633.828 
2098.732 


1.700 
1.700 


120 


0.266 
0.268 


0.29 
0.35 


0.12 
0.13 


























1.5 



150 



10 


2241.796 
2725.402 


1.500 
1.500 


8 


0.266 
0.270 


0.29 
0.70 


0.12 
0.21 


























1.35 



100 



7 


3221.584 
3676.430 


1.350 
1.350 


6 


0.266 
0.269 


0.29 
0.45 


0.12 
0.16 


























1.25 



100 



6 


4352.893 
4680.580 


1.250 
1.250 


4 


0.266 
0.270 


0.29 
0.43 


0.12 
0.15 


























1.1 



50 



3 


5464.714 
5591.467 


1.100 
1.100 


2 


0.266 
0.256 


0.29 
0.30 


0.11 
0.12 


























1 



50 




2 


8540.320 
8788.403 


1.000 
1.000 


2 


0.212 
0.253 


0.29 
0.30 


0.11 
0.12 


























0.9 



50 




2 


13461.100 
13954.450 


0.900 
0.900 


1 


0.206 
0.246 


0.29 
0.30 


0.11 
0.12 


























0.8 



50 



1 


21552.716 
22448.550 


0.800 
0.800 


1 


0.195 
0.235 


0.29 
0.29 


0.11 
0.12 
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The tracks in the HR diagram for the non-rotating models are 
presented in Fig. |4] For the morphology of the tracks, we note 
the following features (going from the massive to the low-mass 
stars): 

- Owing to the effect of mass loss and high luminosity and 
LjM ratio, the MS band widens between 40 and 120 Mq, 
showing a bump around 60 Mq extending to low effective 
temperatures down to log(reif) = 4.1. 

- The maximum luminosity of red supergiants is expected to 
be around log(L/L0) = 5.7. 

- Stars with initial masses between 25 and 40 Mq evolve back 
to the blue after a red supergiant phase. 

- Extended blue loops, crossing the Cepheid instability strip, 
occur for masses between 5 and 9 Mq. 

- The tip of the red giant branch, i.e. the luminosity at which 
the He-flash occurs is around log(L/L0) - 3.4 for stars in 
the mass range between 1.35 and 1.7 Mq. 

In FigurelH the colour code indicates the ratio of the nitrogen 
to the hydrogen number abundances at the surface of the stars, 
on a logarithmic scale (the initial value is log(N/H) H- 12 = 7.78). 
No changes in the surface abundances are expected during the 
MS phase for stars less massive than 60 Mq. Above that mass 
limit, the stellar winds are sufficiently strong to uncover nuclear- 
processed layers already during the MS phase. We note that 
when the hydrogen surface abundance is between 0. 1 and 0.4 (in 
mass fraction), i.e. when the star is in the first part of the WNL 
phase, we have log(N/H) -H 12 ^ 9.2 - 9.8 (green portion of the 
tracks). In those parts, the N/C ratios (in number) are typical of 
the CNO equilibrium and between 70 and 130 (the initial value 
is N/C = 0.25). 

For masses below about 40 Mq, the changes in the surface 
abundances occur after the star has passed through the red su- 
pergiant or red giant stage. The 20 Mq track is on the verge 
of evolving back to the blue. Its maximum N/C ratio is already 
around 70 and its maximum log (N/H)+12 is equal to 9. For 
masses equal to 15 Mq and below, the maximum N/C and log 
(N/H)-i-12 values are below 3 and 8.6, respectively. 

Since the Schaller et al. ( 1992) grids, many physical inputs 
have changed. The initial abundances have been revised, the 
opacities and reaction rates have been updated, the mass loss 
prescriptions are not the same, and the overshoot is different. For 
all these reasons, the comparison of the 92' results with the non- 
rotating models of the present work reveals many differences. 

Figure |5] compares the present non-rotating tracks and our 
previous 92's solar metallicity grid. The most striking difference 
(compare the long-dashed green curve with the dotted red curve) 
is the extension of the MS band which is reduced in the present 
grid as a consequence of the smaller core overshoot (0.10 in the 
present grid instead of 0.20 in the 92's one). The blue loop for the 
present 7 Mq model is more extended and of lower luminosity, 
as expected for a smaller core overshoot and lower initial metal- 
licity. The mass range going through a Cepheid loop has also 
changed: in the 92' grids, the 7 to 12 Mq models pass through 
the Cepheid instability strip, while in the new grid the limits are 
shifted to 5 to 9 Mq. 

As shown in Fig. |6] (dotted red line), the duration of the 
MS is shorter in the new models, by about 15% for the lower 
masses, and 5% for the intermediate masses. Only the most mas- 
sive models (M > 30 Mq) have a longer MS duration by about 
4%, probably owing to the difference in the mass-loss rates. In 
contrast, the duration of the central He-fusion phases is longer 
in the new models, by 33-55% for the intermediate masses and 
10-17% for the massive ones up to 40 Mq. At yet higher masses. 
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Fig. 5. HR diagram of the 3, 7, and 20 Mq models. Comparison 
between the present work rotating models (solid blue line), non- 
rotating models (dashed green line), and the 92' grids models 
(dotted red line). 

Table 3. Coefficients of a piecewise linear fit to calculate the MS 
duration from the mass of the star. 



log(rH [yr]) = A ■ log(M/Mo) + B standard 



mass range 


A 


B 


deviation 


1.25-3.0 


no rot. 


-2.926 


9.892 


0.016 




rot. 


-2.776 


9.938 


0.012 


3-7 


no rot. 


-2.405 


9.641 


0.013 




rot. 


-2.444 


9.761 


0.013 


7- 15 


no rot. 


-1.765 


9.105 


0.010 




rot. 


-1.763 


9.186 


0.015 


15-60 


no rot. 


-0.808 


7.954 


0.025 




rot. 


-0.775 


8.004 


0.022 



the He-burning duration is shorter by 7-26%. These differences 
are mainly due to changes in both the overshooting parameter 
(for masses below 40 Mq) and the mass-loss rates (for masses 
above 40 Mq). 

Figure [T] presents the final masses obta ined for our models 
and the ones given in Schaller et al.l 11991 . We see that for the 
lower mass range, the present models deviate from the slope one 
relation (grey line) at lower mass than the 92's models, typically 
around 15 Mq instead of 25 Mq. This is a consequence of the 
prescription adopted for the mass-loss rate during the red super- 
giant stage (cf. Sect. 12.6.21 ) which leads to stronger mass loss 
during that phase in the present models. 

In the upper mass range (above 35 Mq), the present non- 
rotating models have higher masses at the end of the C-burning 
phase than the masses we obtained in 92. This is a consequence 
of the lower mass-loss rates used here during the Wolf-Rayet 
phases. 
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Fig. 6. Difference in the MS lifetime duration. Comparison be- 
tween the present work rotating and non-rotating models (solid 
blue line), and comparison between the present work's non- 
rotating and the 92' grids models (dotted red line). 



5. Properties of the rotating models 

5.1. Main sequence width 

In Fig. [8] we superimpose the evolutionary tracks of the 
present paper with rotation (solid lines) and without rotation 
(dashed lines) on the data p oints for the intermediat e- mass main- 
sequence stars discussed in lWolff & SimonI (Il997h . We see that 
for masses above about 1 .7 Mq, non-rotating tracks computed 
with an overshoot parameter doyeJHp - 0.10 predict a too nar- 
row MS band, while the tracks computed with an average rota- 
tion and this same overshoot parameter closely fit the observed 
MS band. Thus, in this mass range, the effects of a moderate 
overshoot and rotational mixing closely fit the observations. 

For masses between 1 .25 and 1 .5 Mq, the situation is much 
less clear; the non-rotating MS tracks are in general too narrow 
and the rotating ones too wide. This mass domain just corre- 
sponds to the upper mass range where magnetic braking comes 
into play. We note that with our present prescription for mag- 
netic braking, calibrated for the solar case, the mean velocity 
during the MS phase is between 6 and 10 km s This is smaller 
than the average observed rotations for this mass range, which 
are between 10 and 40 km s"' for field stars and between 20 and 
130 km s"' for cluster stars according to Wolff & Simon ( 1997). 
Hence, it may be that the present models, in this mass range, 
overestimate the magnetic braking, as we will investigate in a 
future paper. Until then, we note that rotating and non-rotating 
models encompass the observed MS widths in this mass range. 

Ide Meulenaer et aO (l2010l) . using asteroseismology, found 
that a Cen A (a G2V star with a mass of 1.105 + 0.007 Mq) has 
observed asteroseismic properties that can be more accurately 
reproduced by models with no convective core confirming the 
choice we made to account for no overshoot in models below 
1.25 Mq. 
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Fig. 7. Final mass vs initial mass for models from 9 to 120 Mq. 
Comparison between the present work's non-rotating models 
(dashed green line) and rotating models (solid blue line), and 
the mass loss rates of the 92' grids models computed with the 
"normal" rates (dotted red line). The grey line corresponds to a 
hypothetical case without mass loss (Mf\„ = Mini). 



With the CoRoT satellite. [Briquet et all ( 1201 lb observed the 
09V star HD 46202 (member of the young open cluster NGC 
2244). The mass of this star is around 20 Mq, its vsini is ~25 
km s and it is on the main sequence (at an age estimated to 
be around 2 Myr). If this star were a truly slow rotator (namely 
if the low V sin / were not due to a sin /-effect and the star were 
throughout its lifetime a slow rotator), then the extension of the 
core would be entirely due to overshooting, with no or an ex- 
tremely small contribution from rotational mixing. Very inter- 
estingly, their best-fit model gives a core overshooting parame- 
ter dover/Hp = 0.10 ± 0.05, quite in agree ment with the over- 
shoot used in the present models. In contrast. [Lovekin & Goupill 
(.2010i) . studying Ophiuchi, a MS ~9.5 Mq star with a u sin ; of 
30 km s'^, obtain an overshooting parameter equal to 0.28+0.05. 
If this star were a truly slow rotator, and thus the extension of the 
core due mainly to the overshooting process, the result obtained 
in this work would be at odds with our choice of overshooting 
parameter. This illustrates that despite many decades of discus- 
sion about the overshooting parameter, its precise value remains 
uncertain. 

5.2. Hertzsprung-Russell diagram and lifetimes 

The tracks in the HR diagram for all the rotating models are 
presented in Fig.|9] with lines of iso-velocities drawn across the 
diagram. Compared to Fig.H] we note the following differences: 

- There is no longer any widening of the MS band occurring 
around the 60 Mq stellar model as in the non-rotating grid. 
This is a consequence of rotational mixing, which prevents 
any significant redward evolution during the MS phase. This 
shows that in this upper mass range, the MS band width is 
very sensitive to rotational mixing. 
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Fig. 9. Same as Fig.|4]but for the rotating models. Lines of iso-velocities are drawn through the diagram. 
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Fig. 8. Evolutionary tracks for non-rotating (dashed lines) and 
rotating (solid lines) models in the colour-magnitude diagram. 
The ZAMS and the minimal MS Teff are indicated as well as the 
initial mass considered. The crosses are observed MS intermedi- 
ate mass stars from Tables 1 and 2 of .Wolff & Simoa(.1997.) . 

- The maximum luminosity of red supergiants is in the range 
log(L/ Lq) = 5 .2 - 5 .4, hence is a factor of two inferior to the 
limit obtained with non-rotating mo dels. This value close ly 
agrees with the upper limit given bv lLevesque et al.l (l2005h . 

- Stars with initial masses between 20 and about 25 Mq evolve 
back to the blue after a red supergiant phase. This range is 
between 25 and 40 Mq in the non-rotating grid. 

- Extended blue loops, crossing the Cepheid instability strip 
occur for masses between 5 and 9 Mq as for the non-rotating 
tracks. For a given mass, the loop however occurs at a higher 
luminosity. The blueward extension is similar or slightly 
shorter. 

- The tip of the red giant branch does not appear to be very 
different in the rotating models. We must stress here that not 
all the models have been computed strictly up to the He- 
flash, preventing us from reaching more definite conclusions. 
The above assessment relies on the case of the 1 .7 Mq model 
that was computed up to that stage. 

- Changes in the surface abundances occur at a much earlier 
evolutionary stage (see Sect. l5.4l below). 

Figure |5] directly compares the rotating and non-rotating 
tracks for the 3,7, and 20 Mq models. We see that the present ro- 
tating tracks end the MS phase at lower eff'ective temperature and 
higher lumino sity than the non-rotat i ng ones as found by pre - 
vious studies (iHeger & Langeill2000l: iMevnet & Maededl200dl) . 
They cross the Hertzsprung gap at a higher luminosity than their 
non-rotating counterparts. This shift in luminosity persists for 
the remainder of the evolution. As a consequence, the Cepheid 
loop of the 7 Mq occurs at a luminosity around 0.15 dex higher 
than in the non-rotating case. We note that this corresponds to 
an evolutionary mass that is about 15% lower for a Cepheid at 
a given luminosity. It has been proposed that mass loss and in- 
ternal mixing might r econcile the evolutionary and puls ational 
masses for Cepheids ( iKelled l2008t iNeilson et all 1201 ll) : rota- 



tion would account for both effects and the masses inferred 
would thus be closer to the pulsational ones. The rotating 20 Mq 
evolves back to the blue in contrast to the non-rotating one, indi- 
cating that rotational mixing, together with enhanced mass loss 
duri ng the red supergiant stage, promote blueward evolution (see 
also lHirschi et al. 2004). 

As shown in Fig. |6](blue solid line), rotation generally en- 
hances the time spent on the MS, by about 25% throughout 
all the masses above 2 Mq. The increase in the MS lifetime 
is caused by rotational mixing, which refuels the core in fresh 
hydrogen. It is striking to see how relatively constant the in- 
crease in the MS lifetime appears to be over this whole mass 
range. This shows that the effect of rotation scales very well with 
the initial Vini/van value. The picture is much less clear for the 
core helium- and carbon-burning lifetimes (not shown in Fig.|6]l. 
Some models have an extended lifetime, others shorter. The net 
result depends on many parameters, such as the choice of diffu- 
sion coefficients (which governs the size of the convective core 
and the mixing from or to the adjacent H- and/or He-burning 
shell), or the mass loss. 

Figure |7] compares the final and initial mass relations. For 
masses lower than 40 Mq, the difference obtained between the 
rotating and non-rotating models is small. Above that mass, the 
rotating models up to 100 Mq lose less mass than the non- 
rotating one. Only in the very upper mass range do the rotating 
models end with lower final masses than the non-rotating ones. 
This is caused by the non-rotating stars spending more time in 
regions of the HR diagram where the mass loss is stronger. By 
remaining in the blue part of the HR diagram, the rotating mod- 
els pass from O-type mass-loss rates directly to WR mass-loss 
rates, without experiencing the bistability jump, and the red su- 
pergiant (RSG) and supra-Eddington mass-loss rates. 

5.3. Evolution of the surface velocities 

Figure [TO]presents the evolution of the equatorial velocity (left) 
and the feq/L'ciit ratio {right), during the MS for all the models. 
Since we computed models with identical initial Vsq/Van ratio, 
the initial equatorial velocity varies with the mass considered: 
higher masses need a larger equatorial velocity to attain a given 
ratio. In Fig. [TOl(/e/r). this is clearly visible: the 1.7 Mq model 
draws the bottom black line, starting its MS evolution with i/gq = 
170 km s"\ while the most massive 120 Mq model draws the 
upper red line on the ZAMS, starting its MS evolution with Weq = 
389kms-i. 

We bring our models onto the ZAMS assuming solid body 
rotation, an assumption that is no longer made after the ZAMS. 
As soon as we allow for differential rotation, a readjustment of 
the Q-profile occurs inside the model, which explains the rapid 
drop in equatorial velocity at the very beginning of the evolu- 
tion (Denissenkov et al. 1999). Once the equilibrium profile is 
reached, the equatorial velocity evolves under the action of the 
internal transport mechanisms (which tend to bring angular mo- 
mentum from the contracting core towards the external layers, 
spinning up the stellar sur facej] and the mass loss (which re- 
moves angular momentum from the surface, spinning down the 
star). 



' Even though the net transport of angular momentum is towards the 
surface, approximating the meridional circulation (an advective pro- 
cess) as a diffusive process would lead to even more transport of an- 
gular momentum from the core to the surface and to less shear mixing 
and thus less mixing of elements like nitrogen to the surface. 
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Fig. 10. Evolution of rotation on the main sequence as a function of the central mass fraction of hydrogen for all the stellar models 
from 1.7 to 120 Mq. Left: Evolution of the equatorial velocity. The curve for the 1.7 Mq is the bottom line at the beginning of the 
evolution (Xc - 0.72), and that for the 120 Mq the upper line (increasing equatorial velocity for increasing mass). Right: Evolution 
of the Veq/i^crit ratio. The curve for the 1.7 Mq is the upper line during the MS, and that for the 120 Mq the bottom line during the 
MS (decreasing velocity ratio for increasing mass). 



In low- and intermediate-mass stars, the transport mecha- 
nisms are quite inefficient, but at the same time, mass loss is 
either absent or very weak (at least during the MS). Stars thus 
evolve with a quas i constant i/gq- This agrees with the results of 
IWolffetal.l(l2007h . who find that the rotation rates for stars in 
the mass range 6-12 Mq do not change by more than 0.1 dex 
over ages between about 1 and 15 Myr. 

With increasing mass, the internal transport of angular mo- 
mentum becomes more efficient, but at the same time the mass 
loss becomes dominant, so the most massive models (M > 
32 Mq) quickly spin down. Therefore we see that the time- 
averaged surface velocity is smaller (see Table 2) for stars in 
the mass range between 32 and 120 Mq (between 100 and 180 
km s"'), than in the mass range from 9 to 25 Mq (~ 200 km s"'). 

It is interesting to compare the mean rotation velocity over 
the whole main sequence vms with the velocity of the stars at 
the end of the MS phase (see Table 2). We see that for stars with 
initial masses equal or superior to 32 Mq, the surface velocity 
at the end of the MS phase is much smaller than vms, reaching 
a minimum of only a few km s"' and reaching at most a value 
of 20 km s"' . For these stars, we observe what might be referred 
to as a very strong mass-loss braking mechanism. The situation 
is qualitatively the same for stars in the mass range between 12 
and 25 Mq, although less spectacular since the mass-loss rates 
are much lower. For the masses below 9 Mq down to 1.7 Mq, 
the situation is different: the velocities at the end of the MS phase 
are nearly equal, and sometimes even slightly superior to vms ■ 
This illustrates, as already indicated above, that in the absence 
of stellar winds, the surface continuously receives angular mo- 
mentum transported by the meridional currents, and to a smaller 
extent by shears. If this were not the case, the surface would slow 
down simply as a result of its inflation and the local conserva- 
tion of the angular momentum. For masses equal to or below 1 .5 
Mq, we see that vms is also very similar to the velocities at the 
end of the MS phase, although much smaller than the initial ve- 



locities. This is the result of the very rapid braking exerted by 
the magnetic fields at the beginning of the evolution. 

After the MS phase, the surface velocities and their evolu- 
tion are depicted by the lines of iso-velocities drawn across the 
diagram in Fig. |9] The surface velocity decreases rapidly dur- 
ing the crossing of the Hertzsprung gap for all masses. For all 
stars between 2 and 12 Mq, the iso-velocity line correspond- 
ing to 10 km s"' is nearly vertical with an abscissa in the range 
log(refF) = 3.65 - 3.7. When a blue loop occurs, the surface ve- 
locity increases, sometimes reaching the critical velocity as in 
the case of the 9 Mq model. In this case, however, the proximity 
to the critical limit occurs over such a short duration that the star 
undergoes very little mechanical mass losses. 

Looking at the velocities obtained at the end of the core He- 
burning phase in Table|2] we see that the stars finishing their lives 
as WR stars have surface velocities between only a few km s"' 
as the 120 Mq model, up to slightly more than 70 km s"' for the 
40 Mq model. This is in good agreement with the observations 
of WR stars bv lChene~& St-Louis (2008), who obtain rotational 
velocities of between 1 and 60 km s ' from periodic large-scale 
wind variability. 

Figure [To] (ng/zf) illustrates that the overall evolution of the 
Veq/vcrit ratio looks diff'erent from that of Ugq, at least for the 
models that do not undergo strong braking by mass loss. This is 
due to the evolution of itself: in the framework of the Roche 
model we use here, we have Ucrit y/M/R. During the MS, M 
decreases with the mass loss, and R inflates, so Vcnt decreases. 
This explains why models of low- and intermediate-mass that 
have a quasi constant or slightly decreasing Ugq have the Ve^/Vcm 
ratio that increases during the MS evolut ion. A more detaile d 
discussion of this subject can be found in Ekstrom et al.l (l2008h . 
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Fig. 12. log(g) - log(reff) diagram {left) and evolution of the abundance ratio N/H (right) for the 7 (blue), 15 (red), and 32 Mq 
(green) rotating models. The track of the non-rotating 15 Mq is indicated with the red d otted line. Observa tional po ints from O- and 
B-type stars are plotted, with empty symbols for luminosity class (LC) V-Ill (triangles: iMorel et al.l l2008: squares: Villamariz et aL 
2002, Villamariz & Herrero 2005; circles: Przybilla et al. 2 010) and fill ed symbols for LC I (triangles: .Takeda & Takada-Hidai 
20001: squares: [Crowtheret al.1120061: hexagons: iSearle et al.ll2^ circles: iPrzvbiUa et al.ll20Tol) . 



5.4. Evolution of tiie surface abundances 

The evolution of the surface abundances depends upon both the 
mixing processes, which modify the composition of the internal 
stellar layers, and the mass loss, which removes the external lay- 
ers and uncovers the deeper regions. Rotation, by driving mixing 
and enhancing mass loss is very efficient in enriching the surface , 
as has be en shown by previous studies ('Heger & Langer"2000^, 
iMevnet & Maeder 2000). The colour code in Fig. |9] illustrates 
this. The variation in the N/H ratio at the surface of various ini- 
tial mass models at two evolutionary stages during the MS phase 
is shown in Fig.[TT] Only models with rotation are shown, since, 
for the mass range considered here, non-rotating models would 
not show any changes in the surface abundances and their track 
would be a horizontal line at an ordinate equal to 0, as can also 
be seen in Fig.|4] 

We see that changes in the surface abundances occur as soon 
as the MS phase for all masses above 1 .25 Mq (see also Table|2l). 
The [N/H] ratio at the middle of the MS is higher by more than 
a factor of two at the surface of all stars with M > 13Mq (see 
Fig.fTTTi. If we consider the end of the MS phase, the minimum 
initial mass for which there is a factor of two enhancement in the 
[N/H] ratio is less than 6 Mq. The variations in the [N/C] ratio 
follow the same trend as the [N/H] ratio. It is simply shifted to 
higher values since the carbon abundance varies faster than the 
hydrogen abundance at the surface of the stars. 

As already discussed in lMaeder et alJ (|2009'), the enrichment 
attained depends on both the initial mass (larger at higher mass) 
and the age (larger at more advanced stages). We see that this 
trend continues down to 1.7 Mq. Below that mass, magnetic 
braking, by enforcing a large gradient in the rotational velocities 
in the outer layers, boosts the mixing of the chemical species (see 
also Meynet et al. 201 1) and produces the peak seen in the low 
mass range in Fig. [TT] As mentioned in Sect. 15.11 the magnetic 
braking might be overestimated in these low-mass models, since 



the predicted enrichment is higher than expected in this mass 
domain. A smaller magnetic braking would probably allow the 
models to closely fit both the observed surface velocities and en- 
richments in this mass range. The present models, however, give 
the consistent solution obtained when a solar magnetic braking 
law is applied to stars with initial masses between 1 and 1.5 Mq. 

Looking at Fig.[TT] we see also that the curves form a kind of 
plateau in the mass range between 15 Mq and 20-25 Mq. This 
plateau is the result of two competing effects: first, as discussed 
in the review by Maeder & Meynet (201 1, in press), for a given 
interior differential rotation, mixing becomes more efficient in 
increasing initial stellar masses; second, when the initial mass in- 
creases, meridional currents also become more rapid in the outer 
layers; they are thus more efficient at eroding the rotational ve- 
locity gradients and make rotational mixing less efficient. Above 
the plateau, the effects of mass loss in uncovering the inner lay- 
ers come into play, and the curves then increase more rapidly. 

The evolution of the surface abundances after the MS phase 
in the rotating models is shown in colour scale in Fig. |9l and 
the N/C and N/O ratios are given in Table |2] When the star first 
crosses the Hertzsprung gap, its surface abundance is the same 
as the one reached at the end of the MS phase, since that phase is 
too short to allow any further changes in the surface abundances. 
As a consequence, slowly or non-rotating MS stars are not ex- 
pected to show any changes in the surface abundances when they 
are on their first track towards the red supergiant stage, while 
those originating from moderately rotating progenitors on the 
MS band, would differ away from the initial abundances. 

During the red supergiant phases, some additional changes 
in the surface abundances occur, resulting from the dredge-up 
and additional rotational mixing. Thus, the stars originating from 
moderately rapid rotators at that stage undergo stronger enrich- 
ments than those originating from either slowly or non-rotating 
stars. Typically, at the end of the core He-burning phase, the N/C 
ratios are of the order of 2 in our non-rotating models with initial 
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Fig. 11. The solid curves show the variation as a function of 
the initial mass of the N/H ratio at the surface at, respectively, 
the end (upper solid blue line labelled by "End of MS") and at 
the middle of the MS phase (lower solid red line labelled by 
"Xc ~ 0.35"). The lower hatched zone corresponds to mean val- 
ues of the N/H excess observed at the surface o f Galactic MS B- 
type stars with initial masses inferior to 20 M(T| ( Gies & Lambert! 
[l99l lKmai3[T99l iMorel et al.ll200i lHunteretal■ll2009^ ■ The 
upper one corresponds to the observed range for maximum val- 
ues. The dashed curves show the variation in the N/C ratio nor- 
maUsed to the value on the ZAMS. 



masses between 2.5 and 15 Mq, while they are between 3 and 
10 in the corresponding rotating models, those of higher masses 
presenting the higher N/C values. The values for the N/O ratio 
are between 0.4 and 0.5 for the non-rotating models and between 
0.6 and 1.3 for the rotating ones. 

As a consequence of the high mass-loss rates experienced by 
the stars between 20 and 32 Mq during the red supergiant phase, 
strong variations in the N/C and N/O ratios are obtained. Values 
as high as those obtained at CNO-equilibrium are reached or in 
contrast, very low values owing to the appearance of He-burning 
products at the surface (because N is destroyed during core He 
burning). 

Above 40 Mq, we see He-burning products at the surface, 
hence the zero values obtained for the N/C and N/O ratios. Only 
the non-rotating 40 Mq model has H-burning products at the 
surface. 

Comparing Fig. |4] with Fig. |9l we see that the portions of 
the tracks for which Alog(N/H) +12 > 10 are more numer- 
ous and longer in the rotating models than in the non-rotating 
ones. This stage corresponds to a transition stage between the 
WN and the WC phases. When rotational mixing is accounted 
fo r, this tran sition is much smoother as already discussed by 
iMevnet & Maedeii (i2003b . and thus this stage is longer. 

In Fig. [I2l the observed positions of the solar metallicity 
stars in the \og{g) versus log(reff), and log(N/H) versus \og(g) 
planes are plotted (see the references for the observations in the 
caption of the figure). To see whether the models provide a good 
match to these observations, one must check that the track for an 
appropriate initial mass closely reproduces both the surface ve- 
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Fig. 13. Evolutionary tracks of the 20 and 25 Mq models in the 
diagram of the log(reff) as a function of the central helium mass 
fraction. In this diagram stars evolve from the left to the right 
during He-burning. 



locity and the surface enrichment at the appropriate evolutionary 
stage. Since for most of the stars only the v sin ; is known (if it 
is known), we do not attempt any comparisons with the surface 
velocities. We merely assume that most of the stars should rotate 
with velocities near the observed averaged one for the B-type 
stars. In that case, the present rotating models should provide a 
reasonably good fit to the averaged surface enrichments. 

The left panel of Fig. [12] indicates that most of the sample 
stars have initial masses between 7 and 32 Mq. In the right panel 
of Fig. [12] we see that globally, the theoretical models closely 
cover the ranges of observed abundances. We note that since the 
sample contains B-type stars, no observed points are plotted be- 
low \og(g) ^ 1 . To facilitate the comparison, we have coloured in 
red the observed stars with the highest N/H values. Interestingly, 
these stars are all evolved stars, and most of them cluster around 
the 32 Mq track in the two panels of Fig.[T2]giving some support 
to the present rotating models. 

For most of the observed points, the error bars in the 
N/H ratios are quite large (of the order of 0.20 dex). Only 
iPrzvbilla et al.l (1201 Oh reported smaller error bars of the order 
of 0.05 dejfl Second, many stars have smaller N/H values than 
the initial value we adopted for our models. This might be be- 
cause in some samples, stars of different initial metallicities, and 
thus different initial N/H ratios, have been observed. This again 
shows that precise and homogeneous samples are needed to pro- 
vide constraints on the stellar models. 

To test our models in a slightly different way, we have su- 
perimposed in Fig. [TT] (see the lower hatched areas) the mean 
values of the N/H enhancements for B-type dwarfs stars with 
masses below 20 Mq as given in Table 2 of Maeder & Meynet 



Note that when the observed points of iPrzvbilla et al.l ( l201(]h are 
plotted in the N/C and N/O diagram they follow the trend expected 
based on the activity of the CNO cycle. Since this trend is not stellar- 
model dependent, but a feature resulting from nuclear physics, it pro- 
vides a strong quality check of the surface abundances determinations. 
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(2011, in pres s). These mean values a re o btained from t he ob- 
servations by iGies & Lamberll (Il992h and iKihanI (11992 ') for a 
sample of 22 stars. The maximum observed values are shown in 
the upper hatched area. We see that the mean observed values 
are closely reproduced by our rotating models at the middle of 
the MS phase, which is quite encouraging. We can also see that 
the maximum values are reproduced well by the situation ob- 
tained in our models at the end of the MS phase. Some extreme 
cases can of course be found for very rapidly rotating models at 
an earlier stage, but the fraction of those stars is probably small 
(except maybe when regions rich in Be-type stars are observed). 

Thus, we conclude from these comparisons that our models 
provide a good fit to the averaged observed surface enrichments. 

5.5. Some consequences for massive star evolution 

We defer a detailed discussion of the Wolf-Rayet populations 
expected from these models to a forthcoming paper (Georgy et 
al., submitted). We discuss here only an interesting consequence 
of the present models for the evolution of stars passing through 
a red supergiant stage. As emphasized above, we implemented 
a new mass-loss prescription for those stars, which significantly 
enhances the mass-loss rate during this phase. We note that this 
enhancement occurs only above a certain mass limit, which is 
between 15 and 20 Mq. This closely corresponds to the interval 
of masses containing the upper initial-mass limit of stars explod - 
ing as type II-P supernovae according to ISmartt et al.l (l2009h . 
These authors indeed find that stars above a mass limit around 
17-18 Mq do not explode as type 11-P supernovae, i.e. do not 
follow the characteristic lightcurve expected in the presence of a 
massive and extended H-rich envelope. Since red supergiants in 
the mass interval between 1 8 and about 30 Mq are observed, one 
might wonder what kind of supernovae stars in this mass range 
would produce. Do they collapse to form a black hole, producing 
faint or even no supernova event? Or do they explode as type II- 
L, or type lb supernova, as a result of the loss of the whole, or at 
least a great fraction, of their H-rich envelope? At the present 
stage, we would clearly favour this second hypothesis . There 
are in deed many observations (Smith et al. 2004; van Loo n et all 
I2OOI indicating that red supergiants (or at least part of them) ex- 
perience stronger mass loss than presently ac counted for in the 
models, which has led recently some authors (lYoon & Cantiellol 
I2OIO ) to explore the physical mechanisms responsible for this. 
In their work, they invoked the possibility that pulsation might 
trigger mass loss. Here, we propose another mechanism, that the 
luminosity reaches the Eddington luminosity in the outer layers. 
As indicated above, this implies that the enhancement will only 
occur above a given mass limit, which seems to be consistent 
with the finding of Smartt et al,, (.2009). 

To predict whether stellar models will explode as either a 
type II-P or a type II-L supernova remains however diflicult. We 
defer this discussion for a forthcoming paper, and only state that 
the present enhanced mass-loss rates favour the explosion as a 
type II-L or even a type lb supernova. 

Stronger mass loss during the red supergiant phase favours a 
bluewards evolution. This is a well known effect obtained many 
times in the past (see for instance Salasnich et al. 1999). Thus, 
it changes the time spent in the blue and the red part of the HR 
diagram du ring the core He-burning phase. In Fig. [13] we see, 
as found bv iMaeder & Mevnetl (12001 ). that the first crossing of 
the Hertzsprung gap occurs more rapidly when rotation is ac- 
counted for. We also see that a large fraction of the end of the 
core He-burning phase (i.e., when Y < 0.30-0.60) occurs in the 
blue part of the HR diagram in the rotating tracks. This is due 



to the combined effects of the enhanced mass-loss rates and ro- 
tational mixing, and has very important consequences for the 
blue to red supergiant ratio (B/R). In the models of Schalle ret al.l 
(I1992I) . the ratio of the blue to the red lifetime (red being defined 
by the time spent with an effective temperature below 3.65) for 
the 20 Mq model is 0.2, while for the present non-rotating and 
rotating model, it is 0.7 and 1.5, respectively. The sam e num- 
bers f or the 25 Mq models are respectively 0.10 ( Sc haller et al.l 
I1992h . 1.6 (present non-rotating), and 5.3 (rotating). Thus in that 
mass range, the B/R ratio is much larger with rotation. Whether 
it is sufficient to explain the B/R ratio observed in clusters re- 
mains to be studied and will be considered in a future paper. We 
note that at the metallici ty of our models, the results found by 
[Eggenbe rger et alj (l2002h are between 1.2 and 4, in close agree- 
ment with our 20-25 Mq rotating models. 

Finally, this enhanced mass loss during the red supergiant 
phase lo wers the minim um mass for single stars to become a 
WR star. lSchaller et alj(ll992) o btained an inferior mass limit of 
32 Mq (see Mevnet et al. 19941) . while the present non-rotating 
and rotating models give 25 and 20 Mq respectively. The conse- 
quences of this result for the WR populations at solar metallicity 
will be discussed in the next paper of this series. 



6. Discussion and conclusion 

We have studied for the first time, the impact of rotation on the 
evolution of single stars covering the mass domain from 0.8 to 
120 Mq in a homogeneous way. The present rotating models 
provide globally a good fit to the observed MS width, to the 
positions of the red giants and red supergiants in the HR dia- 
gram, and to the evolution of the surface velocities and the sur- 
face abundances. It is remarkable that such a good agreement 
can be obtained based on a unique well-chosen value of the ini- 
tial f/wcrit for most of the stellar mass range considered. 

The rotating models provide a more realistic view than non- 
rotating models in the sense that only models including rota- 
tional mixing are able to account for all the above features at 
the same time. Using a larger overshoot in non-rotating models, 
for instance, could reproduce the observed MS width but not the 
observed changes in the surface abundances and obviously not 
the surface velocities. All the data as well as the corresponding 
isochrones are made available through the wet0. 

In addition to an extended grid of models with rotation, this 
work also illustrates the consequences of our revised mass-loss 
rate for the red supergiant phase. With this prescription, the stars 
between 20 and 32 Mq lose very large amounts of mass, en- 
abling some of them to evolve back to the blue part of the HR 
diagram, to change significantly the times spent in the blue and 
the red part of the HR diagram during the core He-burning phase, 
to lower the initial mass of stars that become WR stars and to 
change the lightcurve at the time of the supernova event. 

In the future, we will provide grids for other initial metal- 
licities. A very fine grid in mass for low-mass stars (without ro- 
tation) has been computed by Mowlavi et al. (submitted). This 
grid, computed with the same physics as here, allows very accu- 
rate interpolations for both age and mass determinations. 

Acknowledgements. The authors express their gratitude to Dr J. W. Ferguson, 
who has computed on request the molecular opacities for the peculiar mixture 
used in this paper They also thank Claus Leitherer for his careful checking of the 
electronic tiles and preliminary work with the new models. R. H. acknowledges 



' |http : //obswww . unige . ch/Recherche/e vol/-Da tabase- 1 
or through the CDS database at 

http : //vizier .u- strasbg . fr/viz-bin/VizieR-2 | 



Ekstrom et al.: Grids of stellar models with rotation 



19 



support from the World Premier International Research Center Initiative (WPI 
Initiative), MEXT, Japan. T.D. acknowledges support from the Swiss National 
Science Foundation (SNSF), from ESF-EuroGENESIS, and from the French 
Programme National de Physique Stellaire (PNPS) of CNRS/INSU. 



References 

Alpar, M. A., Ankay, A., & Yazgan, E. 2001, ApJ, 557, L61 

Angulo, C, Amould, M., Rayet, M., et al. 1999, Nucl. Phys. A, 656, 3 

Angulo, C. & Descouvemont, P 2001, Nucl. Phys. A, 690, 755 

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in ASPC, Vol. 336, Cosmic 

Abundances as Records of Stellar Evolution and Nucleosynthesis, ed. T. G. 

Barnes, III & F N. Bash (San Francisco: ASP), 25 
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Bohm-Vitense, E. 1981, ARA&A, 19, 295 
Bouret, J.-C, Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 
Briquet, M., Aerts, C, Baglin, A., et al. 2011, A&A, 527, A112 
Brown, T. M., Christensen-Dalsgaard, J., Dziembowski, W. A., et al. 1989, ApJ, 

343, 526 

Cacciari, C. 2009, Mem. Soc. Astron. Italiana, 80, 97 
Chaboyer, B. & Zahn, J.-P 1992, A&A, 253, 173 

Chene, A.-N. & St-Louis, N. 2008, in RevMexAA Conference Series, Vol. 33, 
Massive Stars: Fundamental Parameters and Circumstellar Interactions 
(Revista Mexicana de Astronomia y Astrofisica), 32 

Chiappini, C, Hirschi, R., Meynet, G., et al. 2006, A&A, 449, L27 

Couvidat, S., Garcia, R. A., Turck-Chieze, S., et al. 2003, ApJ, 597, L77 

Crowther, P A. 2000, A&A, 356, 191 

Crowther, P A., Lennon, D. J., & Walborn, N. R. 2006, A&A, 446, 279 

Cunha, K., Hubeny, I., & Lanz, T. 2006, ApJ, 647, L143 

Cunha, M. S., Aerts, C, Christensen-Dalsgaard, J., et al. 2007, A&A Rev., 14, 

217 

de Jager, C, Nieuwenhuijzen, H., & van der Hucht, K. A. 1988, A&AS, 72, 259 
de Meulenaer, P, Carrier, F., MigUo, A., et al. 2010, A&A, 523, A54 
Decressin, T., Meynet, G., Charbonnel, C, Prantzos, N., & Ekstrom, S. 2007, 
A&A, 464, 1029 

Denissenkov, P A., Ivanova, N. S., & Weiss, A. 1999, A&A, 341, 181 
Denissenkov, P. A. & Pinsonneault, M. 2007, ApJ, 655, 1157 
Dufton, P L., Smartt, S. J., Lee, J. K., et al. 2006, A&A, 457, 265 
Eddington, A. S. 1925, The Observatory, 48, 73 
Eggenberger, P, Maeder, A., & Meynet, G. 2005, A&A, 440, L9 
Eggenberger, P, Meynet, G., & Maeder, A. 2002, A&A, 386, 576 
Eggenberger, P., Meynet, G., Maeder, A., et al. 2008, Ap&SS, 316, 43 
Ekstrom, S., Coc, A., Descouvemont, P, et al. 2010, A&A, 514, A62 
Ekstrom, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 
Ferguson, J. W., Alexander, D. R., Allard, F, et al. 2005, ApJ, 623, 585 
Flower, R J. 1977, A&A, 54, 31 

Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 
Fynbo, H. O. U., Diget, C. A., Bergmann, U. C, et al. 2005, Nature, 433, 136 
Georgy, C. 2010, PhD thesis, Geneva Observatory (Geneva University) 
Georgy, C, Meynet, G., & Maeder, A. 201 1, A&A, 527, A52 
Gies, D. R. & Lambert, D. L. 1992, ApJ, 387, 673 
Grafener, G. & Hamann, W.-R. 2008, A&A, 482, 945 

Hale, S. E., Champagne, A. E., Iliadis, C, et al. 2002, Phys. Rev. C, 65, 015801 

Heger, A. & Langer, N. 2000, ApJ, 544, 1016 

Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 

Herwig, F, Austin, S. M., & Lattanzio, J. C. 2006, Phys. Rev. C, 73, 025802 

Hirschi, R., Meynet, G., & Maeder, A. 2004, A&A, 425, 649 

Huang, W. & Gies, D. R. 2006, ApJ, 648, 580 

Huang, W., Gies, D. R., & McSwain, M. V. 2010, ApJ, 722, 605 

Hunter, I., Brott, I., Langer, N., et al. 2009, A&A, 496, 841 

Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 

lUadis, C, D' Auria, J. M., Starrfield, S., Thompson, W. J., & Wiescher, M. 2001, 

ApJS, 134, 151 
lUarionov, A. F & Sunyaev, R. A. 1975, A&A, 39, 185 
Imbriani, G., Costantini, H., Formicola, A., et al. 2004, A&A, 420, 625 
Imbriani, G., Limongi, M., Gialanella, L., et al. 2001, ApJ, 558, 903 
Itoh, N., Adachi, T., Nakagawa, M., Kohyama, Y., & Munakata, H. 1989, ApJ, 

339, 354 

Itoh, N., Hayashi, H., Nishikawa, A., & Kohyama, Y. 1996, ApJS, 102, 411 
Jaeger, M., Kunz, R., Mayer, A., et al. 2001, Phys. Rev. Lett., 87, 202501 
Karakas, A. I., Lugaro, M. A., Wiescher, M., Gorres, J., & Ugalde, C. 2006, ApJ, 
643,471 

Kawaler, S. D. 2009, in ASPC, Vol. 416, Solar-Stellar Dynamos as Revealed 
by Helio- and Asteroseismology, ed. M. Dikpati, T. Arentoft, I. Gonzalez 
Hernandez, C. Lindsey, & F. Hill (San Francisco: Astronomical Society of 
the Pacific), p. 385 

Keller, S. C. 2008, ApJ, 677, 483 



Kilian, J. 1992, A&A, 262, 171 

Kippenhahn, R., Weigert, A., & Hofmeister, E. 1967, Methods of Computational 
Physics, Vol. 7, Methods for Calculating Stellar Evolution (New York: 

Academic Press) 

Kosovichev, A. G., Schou, J., Scherrer, P H., et al. 1997, Sol. Phys., 170, 43 
Krishnamurthi, A., Pinsonneauh, M. H., Barnes, S., & Sofia, S. 1997, ApJ, 480, 

303 

Kunz, R., Fey, M., Jaeger, M., et al. 2002, ApJ, 567, 643 
Langer, N. & Maeder, A. 1995, A&A, 295, 685 

Lemut, A., Bemmerer, D., Confortola, F., et al. 2006, Phys. Lett. B, 634, 483 
Levesque, E. M., Massey, P, Olsen, K. A. G., et al. 2005, ApJ, 628, 973 
Lodders, K. 2003, ApJ, 591, 1220 
Lovekin, C. C. & Goupil, M.-J. 2010, A&A, 515, A58 

Maeder, A. 1987, A&A, 173, 247 
Maeder, A. 1997, A&A, 321, 134 
Maeder, A. 1999, A&A, 347, 185 

Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Berlin 

Heidelberg: Springer Verlag) 
Maeder, A. & Meynet, G. 1987, A&A, 182, 243 
Maeder, A. & Meynet, G. 1989, A&A, 210, 155 
Maeder, A. & Meynet, G. 2000, A&A, 361, 159 
Maeder, A. & Meynet, G. 2001, A&A, 373, 555 
Maeder, A. & Meynet, G. 2005, A&A, 440, 1041 
Maeder, A. & Meynet, G. 2010, New A Rev., 54, 32 

Maeder, A., Meynet, G., Ekstrom, S., & Georgy, C. 2009, Communications in 

Asteroseismology, 158, 72 
Maeder, A. & Zahn, J.-P 1998, A&A, 334, 1000 

Malagnini, M. L., Morossi, C, Rossi, L., & Kuracz, R. L. 1986, A&A, 162, 140 

Meylan, G. & Maeder, A. 1982, A&A, 108, 148 

Meynet, G., Eggenberger, P, & Maeder, A. 2011, A&A, 525, LI I 

Meynet, G. & Maeder, A. 1997, A&A, 321, 465 

Meynet, G. & Maeder, A. 2000, A&A, 361, 101 

Meynet, G. & Maeder, A. 2003, A&A, 404, 975 

Meynet, G., Maeder, A., Schaller, G., Schaerer, D., & Charbonnel, C. 1994, 

A&AS, 103, 97 
Morel, T, Hubrig, S., & Briquet, M. 2008, A&A, 481, 453 
Muijres, L. E., de Koter, A., Vink, J. S., et al. 201 1, A&A, 526, A32 
Mukhamedzhanov, A. M., Bern, P, Brown, B. A., et al. 2003, Phys. Rev. C, 67, 

065804 

Neilson, H. R., Cantiello, M., & Langer, N. 201 1, A&A, 529, L9 

Nomoto, K., Tominaga, N., Umeda, H., et al. 2005, in ASPC Series, Vol. 332, 

The Fate of the Most Massive Stars, ed. R. Humphreys & K. Stanek (San 

Francisco: Astronomical Society of the Pacific), 374 
Nugis, T. & Lamers, H. J. G. L. M. 2000, A&A, 360, 227 
Palacios, A., Talon, S., Charbonnel, C, & Forestmi, M. 2003, A&A, 399, 603 
Przybilla, N., Firnstein, M., Nieva, M. F, Meynet, G., & Maeder, A. 2010, A&A, 

517, A38 

Reimers, D. 1975, Memoires of the Societe Royale des Sciences de Liege, 8, 369 

Reimers, D. 1977, A&A, 61, 217 

Rogers, F J. & Nayfonov, A. 2002, ApJ, 576, 1064 

Salasnich, B., Bressan, A., & Chiosi, C. 1999, A&A, 342, 131 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 

Schatzman, E. 1993, A&A, 279, 431 

Schmidt- Kaler, T. 1982, in The Landolt-Bornstein Database - Numerical Data 
and Functional Relationships in Science and Technology, Vol. 2b, Stars and 
Star Clusters, ed. K. Schaifers & H. H. Voigt (Berlin/Heidelberg: Springer- 
Verlag), 451 

Searle, S. C, Prinja, R. K., Massa, D., & Ryans, R. 2008, A&A, 481, 777 
Sekiguchi, M. & Fukugita, M. 2000, AJ, 120, 1072 

Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 
395, 1409 

Smith, N., Vink, J. S., & de Koter, A. 2004, ApJ, 615, 475 
Stothers, R. & Chin, C.-W. 1973, ApJ, 179, 555 

Suarez, J. C, Garrido, R., & Goupil, M. J. 2007, Communications in 

Asteroseismology, 150,211 
Suijs, M. P L., Langer, N., Poelarends, A.-J., et al. 2008, A&A, 481, L87 
Sweet, P A. 1950, MNRAS, 110, 548 

Sylvester, R. J., Skinner, C. J., & Barlow, M. J. 1998, MNRAS, 301, 1083 
Takeda, Y. & Takada-Hidai, M. 2000, PAS J, 52, 113 
Talon, S. & Charbonnel, C. 2005, A&A, 440, 981 
Talon, S., Kumar, P, & Zahn, J.-R 2002, ApJ, 574, L175 
Talon, S. & Zahn, J.-P 1997, A&A, 317, 749 

Talon, S., Zahn, J.-P, Maeder, A., & Meynet, G. 1997, A&A, 322, 209 

Tammann, G. A., Sandage, A., & Reindl, B. 2003, A&A, 404, 423 

The, L.-S., El Eid, M. F, & Meyer, B. S. 2007, ApJ, 655, 1058 

Thompson, T. A., Chang, R, & Quataert, E. 2004, ApJ, 61 1, 380 

Thoul, A. 2009, Communications in Asteroseismology, 159, 35 

Townsend, R. H. D., Oksala, M. E., Cohen, D. H., Owocki, S. P, & ud-Doula, 



20 



Ekstrom et al.: Grids of stellar models with rotation 



A. 2010, ApJ, 714, L318 
Tur, C, Heger, A., & Austin, S. M. 2009, ApJ, 702, 1068 
Tur, C, Heger, A., & Austin, S. M. 2010, ApJ, 718, 357 
Turck-Chieze, S., Piau, L., & Couvidat, S. 2011, ApJ, 731, L29 
van Loon, J. T., Cioni, M.-R. L., Zijlstra, A. A., & Loup, C. 2005, A&A, 438, 

273 

van Loon, J. T., Groenewegen, M. A. T., de Koter, A., et al. 1999, A&A, 351, 

559 

Vazquez, G. A., Leitherer, C, Schaerer, D., Meynet, G., & Maeder, A. 2007, 

ApJ, 663, 995 
Villamariz, M. R. & Herrero, A. 2005, A&A, 442, 263 

Villamariz, M. R., Herrero, A., Becker, S. R., & Butler, K. 2002, A&A, 388, 940 
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 
von Zeipel, H. 1924, MNRAS, 84, 665 

Walder, R., Folini, D., & Meynet, G. 2011, Space Sci. Rev., 57 

Weaver, T. A. & Woosley, S. E. 1993, Phys. Rep., 227, 65 

Weiss, A., Serenelli, A., Kitsikis, A., Schlattl, H., & Christensen-Dalsgaard, J. 

2005, A&A, 441, 1129 
Wolfi; S. & Simon, T. 1997, PASP, 109, 759 

Wolff, S. C, Strom, S. E., Dror, D., & Venn, K. 2007, AJ, 133, 1092 
Woosley, S. E. & Heger, A. 2004, in lAU Symposium, Vol. 215, Stellar Rotation, 

ed. A. Maeder & P. Eenens (San Francisco: Astronomical Society of the 

Pacific), 601 

Yoon, S., Langer, N., & Norman, C. 2006, A&A, 460, 199 
Yoon, S.-C. & Cantiello, M. 2010, ApJ, 717, L62 
Zahn, L-P 1992, A&A, 265, 115 

Zahn, J.-P, Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 
Zahn, J.-P, Talon, S., & Matias, J. 1997, A&A, 322, 320 



